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ABSTRACT 



We consider the possibility that a mass of ~ 10~ 5 — 10~ 3 M Q flows back 
from the dense shell of planetary nebulae and is accreted by the central star 
during the planetary nebula phase. This backflowing mass is expected to have a 
significant specific angular momentum even in (rare) spherical planetary nebulae, 
such that a transient accretion disk might be formed. This mass might influence 
the occurrence and properties of a very late thermal pulse (VLTP), and might 
even trigger it. For example, the rapidly rotating outer layer, and the disk if 
still exist, might lead to axisymmetrical mass ejection by the VLTP. Unstable 
burning of accreted hydrogen might result in a mild flash of the hydrogen shell, 
also accompanied by axisymmetrical ejection. 



1. INTRODUCTION 

A planetary nebula (PN) is an expanding ionized circumstellar cloud that was ejected 
during the asymptotic giant branch (AGB) phase of the stellar progenitor. The nebula is 
ionized by the hot central star. After several thousands years the nuclear reactions in the 
central star cease, when more or less that star starts its cooling track toward the the WD 
track (e.g., Kovetz & Harpaz 1981; Schonberner 1981, 1983). The PN shell is very dense, 
while the region between the central star and the shell has a very low density. One of the 
questions at the early history of PN modelling was why part of the shell mass does not 
flow back to toward the central star? This was answered with the suggestion (Kwok et al. 
1978) and discovery (Heap 1979; Cerruti-Sola & Perinotto 1985) of a fast and tenuous wind 
blown by the central stars of PNs. The fast wind exerts pressure on the PN shell, mainly by 
forming a hot low density bubble, accelerating the shell outward, and preventing mass from 
flowing back. 
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However, under some circumstances small amounts of mass might flow back toward the 
central star. Accretion of backflowing material during the post-AGB phase was mentioned 
before by Mathis & Lamers (1992), Bujarrabal et al. (1998), and Zijlstra et al. (2001), 
and studied by Soker (2001). Most of these studies were interested in substantial mass, 
~ 0.001 — 0.1M o , backflowing from the nebular shell and accreted by the central star. A 
much smaller mass, ~ 1O _6 M , is required in the scenario, involving segregation of gas and 
dust in the outflow, invoked by Mathis & Lamers (1992) to explain chemical abundance 
peculiarities observed in some post-AGB stars. 

The backflowing mass might influence the evolution by forming accretion disk and jets, 
and by supplying fresh material to the core and by that extending the post-AGB (pre-PN) 
phase. Soker (2001) showed that such a process requires that some of the material ejected 
during the AGB phase has a very low outflow speed < 1 km s _1 , such that it can slow down 
and fall within a few thousand years. According to Soker (2001), the back flowing material 
falls because of the gravitational attraction of the central star (or binary system), and falls 
from a distance of ~ 400 AU. A binary interaction can lead to such an outflow. 

In the present paper we consider a different kind of backflow, where gravity is not 
important. We consider a relatively small mass, ~ 10~ 5 — 10~ 3 M Q , that falls back during 
the PN phase rather than the pre-PN phase, and is pushed back to the center by the thermal 
pressure of the ionized shell rather than by gravity. This is discussed in section |2j 

The small amount of mass that falls on the central star at the PN phase might influence 
a very late thermal pulse (VLTP), or even trigger it. Hajduk et al. (2007) considered the 
possibility that mass accretion can induce a VLTP. They discuss this for the 'old nova' CK 
Vul, where the mass is assumed to come from a companion. We instead, discuss the case 
where the accreted mass originates in the nebula (Sect. [2]). 

A VLTP is a helium shell flash occurring so late in the post-AGB evolution that the 
hydrogen burning shell is already extinct (e.g., Fujimoto 1977; Iben et al. 1983; Herwig et al. 
1999). This allows the H-rich matter to be ingested directly into the convective flash zone 
and leads to a rapid, looped evolution on H-R diagram (Herwig et al. 1999). Other "late" 
thermal pulse scenarios, occurring earlier in the course of the post-AGB evolution, before 
complete exhaustion of the H-burning shell, are also considered in Sect. [3j but it appears 
that VLTP pertains the most to our fall-back scenario. The possible influence of the accreted 
mass on the VLTP is discussed in section |3j A short summary follows in section HI 
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2. BACKFLOW DURING THE PLANETARY NEBULA PHASE 

2.1. Formation of backfiowing blobs 

We envision the following process that might lead to the formation of blobs that flow 
toward the central star. 

The initial interaction of the fast wind blown by the post-AGB star with the dense 
shell blown during the end of the AGB phase is momentum-driven, since the shocked fast 
wind material cools very fast and the acceleration of the shell is determined by momentum 
conservation along each direction. When the fast wind speed is high enough the post shock 
radiative cooling time is longer than the flow time and a hot bubble is formed. Radiative 
cooling rate is small, and most of the energy is retained in the flow. The shell is in the 
energy-driven phase; see early works by Volk & Kwok (1985) and Kahn & Breitschwerdt 
(1990). The transition from momentum-driven to energy-driven shell occurs when the fast 
wind speed is Vf ~ 160 km s" 1 (Kahn & Breitschwerdt 1990). The interaction is prone to 
instabilities: the thin shell instability during the momentum-driven phase (Dwarkadas & 
Balick 1998), and the Rayleigh- Taylor instabilities during the energy-driven phase (Kahn & 
Breitschwerdt 1990; Dwarkadas & Balick 1998). 

The energy-driven phase is likely to start before the ionization phase starts, or at about 
the same time. To show this we notice the following relation between the escape speed from 
the central star t> esc , the central star effective temperature T*, its mass M*, and its luminosity 

v- = 170 ( T ' ) (-*-) {^-Y" km s-V (1) 

vi.5xio 4 Ky v°- 6M ©/ V 3oo °W 

Assuming that the wind speed is Vf ~ t> csc , we find that the transition to the energy-driven 
shell occurring at Vf ~ 160 — 200 km s _1 (Kahn & Breitschwerdt 1990), happens at about the 
same time that ionization starts to be significant. This consideration is in accord with the 
central star model of Perinotto et al. (2004), who used a star with luminosity of L* = 6300L Q , 
and where the wind speed of Vf ~ 200 km s -1 is reached when the central star temperature 
is T* = 1.6 x 10 4 K. 

The order of events implies that the hot bubble develops before the dense shell of AGB 
wind is fully ionized, and with it the instabilities (Kahn & Breitschwerdt 1990; Dwarkadas 
& Balick 1998). Before ionization starts the pressure in the dense shell is of the order of 
the pressure in the hot bubble. When ionization becomes important it rapidly propagates 
through the shell, and heats it to a temperature of ~ 10 4 K (Perinotto et al. 2004). The shell 
pressure increases because both temperature and the number density of particles increase. 
The temperature increases from several x 100 K to ~ 10 4 K, so we can take a pressure jump 
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of ~ 100. The shell starts to expand into the hot bubble. Because the initial pressure of the 
ionized shell is much larger than that of the counter pressure of the hot bubble, the front of 
the shell expands at a velocity v exp relative to the shell given by (Landau & Lifshits 1987, 
§99) v cxp ~ 2co/(7 — 1)5 where 7 is the adiabatic exponent and Co — 10 km s -1 is the sound 
speed in the ionized shell. For 7 = 5/3 and 7 = 1.1 one finds v cxp = 3co and v cxp = 20co, 
respectively. 

The ionizing radiation will keep the inflowing mass at a temperature of ~ 10 4 K; a 
value of 7 close to unity is appropriate for the flow. The front of the flow does not contain 
much mass, and therefore the relevant flow speed is a lower velocity, say (0.1 — 0.5) x t> cxp = 
(0.2 — 1) x Co/ (7 — 1) — 20 — 100 km s _1 , where we took 7 ~ 1.1. However, the acceleration 
will be less efficient than the one dimensional flow derived in Landau &Lifshits (1987, §99), 
because any blob formed by instability will expand to the sides (transversely) as well. 

In any case, we conclude that dense blobs formed by instability will start to flow inward 
at a speed of ~ 10 — 50 km s -1 . This speed is larger than the expansion velocity of the 
shell ~ 10 — 20 km s -1 , and there will be a backflow at speed f bac k ~ 3 — 30 km s _1 made 
up of blobs. In the ID calculation of Perinotto et al. (2004) the initial expansion velocity 
of the shell is 10 km s -1 . After ionization starts a zone to the main shell is formed with 
an expansion velocity of only ~ 4 km s -1 . Namely, the inner zone of the main shell was 
substantially slowed down. The pressure of the hot bubble prevents further acceleration 
inward, and later accelerates the entire shell outward. However, if a blob is disconnected 
from the main shell, it will have the hot bubble material all around it. The pressure of the 
hot bubble is more or less isotropic and it will not slow the blob down. This effect cannot 
be studied with a ID code. 

The end product of the initial evolution is a large number of dense (relative to the hot 
bubble) blobs at a temperature of ~ 10 4 K that are immersed in the hot low density gas 
of the hot bubble. Many of these will be flowing inward due to inward acceleration by the 
pressure of the ionized shell during the early ionization stage. If a blob does not slow down 
much, its backflow time from an initial radius r is 

-^-000 ( I5 ^_)( I ^ I )- y , (2 ) 

This is about the time required for the mass to be accreted during the PN phase, but before 
a VLTP occurs. 
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2.2. The properties of the blobs 

For the blob not to slow much it should encounter a mass M e not much larger than its 
initial mass Mb = 47rafpb/3, where a& is the radius of the blob and pb is its density. For 
simplicity we assume the blob is in the shape of a sphere. After flowing a distance ~ r the 
blob encounters a mass 

M e ~ ira 2 b r p h , (3) 

and the condition reads 

l<^40o(W-*-Y (4) 

The density ratio of 300 comes from pressure equilibrium between the hot bubble and the 
blob, and a temperature ratio of ~ 3 x 10 6 K/10 4 K = 300. The blobs' radius is therefore 
o-b > 0.0025r The corresponding mass of the blob is 

Mb ~ M ncb f^X = 10- 7 ( ( M (5) 

where M ne b is the mass in the nebula, and in the first equality we assume that the blob 
density is about equal to the density in the nebula. 

To survive, the blob must be thermally isolated from the surrounding (otherwise it will 
be evaporated). We assume that tangled magnetic fields within the blob isolate the surviving 
blobs. We cannot estimate the number of surviving blobs, but the requirements of (i) an 
initial perturbation that will form a blob with negative radial velocity of vt, ~ 10 km s -1 , 
(ii) a large enough blob to maintain its speed, and (Hi) suppression of heat conduction by 
magnetic fields, suggest that their number will be limited. We scale the number of blobs by 
N b = 1000. 



2.3. Constraints on the fast wind evolution 

As long as the dense blob flows inward inside the hot bubble it does not feel the ram 
pressure of the fast wind. When it gets out of the hot bubble through the inner boundary 
(reverse shock), the fast wind hits the blob directly, exerting a force of F ram = pjvjiral, 
where p/ is the density of the fast wind. The deceleration is a ram = F ram /M b . Using this 
value we find the distance over which the blob will stop, r stop — v b 2 /2a ram , to be 



- 6- 



While inside the hot bubble the density in the blob is determined by pressure equilibrium 
with the hot bubble. Close to the inner boundary of the hot bubble, r bubble , the pressure is 
about equal to the ram pressure of the fast wind there, and so is the blob pressure. This 
gives pb — Pf(vf/c ) 2 , where c — 10 km s -1 is the sound speed in the blob. In our scenario 
Co — Vb, and we find from equation ([6]) the condition for the blob to reach the center to be 
a>b ^ ^bubble- Namely, the blob must be large when it gets out of the hot bubble; we'll come 
back to this point later. 

Equating the ram pressure to the pressure in the blob p/t>? ~ p^c 2 ,, and substituting for 
the fast wind density pf = M//(47rr 2 t>/), and for the density in the blob Mj,/ (47^/3), we 
can isolate a^- Inserting this expression for the blob radius a& in the condition > r bubble 
we derive the condition 



-1 



where in the second equality we have substituted Co = 10 km s 

When the fast wind properties are constant with time, the inner shock (the inner bound- 
ary of the hot bubble r bubb i e ) moves outward as the nebula expands. We require, therefore, 
that within a short time the fast wind momentum discharge, defined as MfVf, will decrease 
by a large factor. In that case r bubb i e oc (M/t>/) 1//2 . In the model of Perinotto et al. (2004), 
Mf drops from ~ 1O~ 8 M to ~ 10~ 10 Af o vr -1 within ~ 1000 yr at an age of t ~ 7000 yr. 
The fast wind speed then is ~ 10 4 km s -1 . Taking values for the inner hot bubble boundary 
from their model, as well as other parameters, we crudely scale the parameters as 

ti — ~r 7 ) ( — -V 2 cm. (8) 

10- n M Q yr- 1 J V4000 km s" 1 / V ; 

Substituting this crude relation in equation ([7]) we find for our demand for mass accretion 
to occur 



Mj \ ( v f \ ( M h x 2/3 



lO- n M yr- 1 / V4000 km s" 1 / ~ V. 10_7M g 



(9) 



The condition at, > r bubb i e for the mass to reach the center as derived above, implies that 
the backflowing gas covers a large solid angle. Let the solid angle covered by the backflowing 
gas be 47T/3. The ratio of the force due to the ram pressure, Fj = f3MfVf of the fast wind to 
the gravitational force of the central star Fq = GM^Mb/r 2 is 

El = ( r \ 2 ( Mi ) ( v _i \ ( M * \ 1 f^L \ 1 (JL 

F g V5 x 10 14 cm/ I 10- n M o yr- 1 I V4000 km s- 1 / V 1O_7M 0/ V°- 6M 0/ V - 25 

(10) 
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From the derivations above we learn the following. Equations (Q shows that for the 
backflowing blob to reach the center the fast wind should be extremely weak, e.g., the mass 
loss rate is an order of magnitude below that in the model of Perinotto et al. (2004) at the 
age of 8000 years. Considering the many unknown in the evolution of the central stars of 
PNs, we consider this constraint on the fast wind not unreasonable. Namely, we conjecture 
that it applies in many cases at the age of several thousand years. From equation (fit)]) it turns 
out that gravity becomes important when the backflowing warm gas reaches a distance of 
~ 100 AU under these conditions. The constraint on the blob mass is Mb > 10 -7 — 10 _6 M Q . 



2.4. Angular momentum 

2.4-1 ■ Single central star 

In section 2.1 we considered the radial acceleration by the nebular thermal pressure of 
small amount of mass toward the center. Let us assume that because of stochastic motion 
in the unstable layer the pressure gradient exerts a small azimuthal (tangential) acceleration 
in addition to the radial acceleration inward. Let the typical azimuthal velocity be 5vb ~ 
0.01(5/0.001) km s _1 , such that the typical angular momentum is jb = 5v},r . The combined 
specific angular momentum of Nb blobs with randomly oriented angular momentum is j ~ 

— 1/2 

N b jb- This angular momentum corresponds to a Keplerian disk at a radius around the 
central star of mass M* given by 

This radius is larger than the radius of the central star during the PN phase, e.g., the VLTP 
occurs when the radius of the central star might be as small as ~ 10 9 cm (e.g., Iben et 
al. 1983). We conclude that an accretion disk might be formed around the central star. 
The accretion is likely to occur over hundreds to thousands of years, with an average mass 
accretion rate of M acc ~ 10~ 9 — 10 _6 Af Q yr _1 . Jets might be launched. In any case, if the 
disk exists during the VLTP, the mass ejection, perhaps resembling a nova eruption, will be 
axisymmetrical. If the PN is not spherical, then the nebula itself might have some amount 
of initial angular momentum, making an accretion disk formation more likely. 



2-4-2. Interaction with a wide companion 



An interaction with a wide companion is another process by which the backflowing 
material can acquire angular momentum. We will take the companion to be of low mass, 
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e.g., a low mass main sequence star, a brown dwarf, or a massive planet. This process is more 
general, and can be relevant to inflow from smaller radii, and at earlier phases of evolution. 
For example, during the final AGB and early post-AGB phases if mass falls back from the 
'effervescent zone' (Soker 2008). The effervescent zone is an extended zone conjectured to 
exist above evolved AGB stars (those with high mass loss rates) extending to ~ 100—300 AU. 
In addition to the escaping wind, in this zone there are parcels of gas that do not reach the 
escape velocity. These parcels of dense gas rise slowly and then fall back. 

Let the companion orbit the central star at radius a. Gas falling to the center from 
radius tq > a will have an infall velocity of 



2GM* ( _ a_ 
r 



1/2 

(12) 



The orbital velocity of the companion is v 2 = (GM^/a) 1 ^ 2 . The relative velocity of the 
inflowing gas and the companion is given by v 2 eX = v\ + vf n ~ 2GM*/a, where we assume 
that r > 1.5a. The accretion radius of the companion is the Bondi-Hoyle-Lyttleton one 

2GM 2 M 2 

where M 2 is the mass of the companion. 

The mass with an impact parameter of b < R a2 is accreted by the companion. Mass at 
larger distances applies drag force on the secondary, with a magnitude of (Alexander et al. 
1976) 

F = 7rR 2 a2 pv 2 el HR max /R a2 ). (14) 

The same force is applied by the star on the gas. The azimuthal component of this force 
is F d = Fv 2 /v Te \. We also take the maximum radius of influence to be -R max — which by 
equation ([13]) gives ln(R max /R a2 ) ~ m(M*/M 2 ) ~ 1, for M 2 ~ 0.1M* ~ O.1M . 

The angular momentum imparted to the inflowing gas by this gravitational interaction 

is 

— = F e a ~ TcR 2 a2 pv re iv 2 a. (15) 

The total inflow rate through radius a is M- m = Aiia 2 pVi n . The specific angular momentum 
of the inflowing gas is j = (dJ / dt) / M- m . Using our approximations, the value of j is 

Jin " {w) V2a - (16) 

To examine the possibility of disk formation around the central star, j in should be compared 
with the specific angular momentum of a Keplerian circular orbit on the star equator = 
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(GM^i?*) 1 / 2 . We find, after using the expression for v 2 , 

j* \Mj \rJ \0.1MJ VlOOAU/ \IRqJ v ; 

where in the second equality we have substituted values appropriate for a central star of a 
PN in our studied case of a wide low mass companion. 

We see that a binary companion can spin-up the inflowing gas to a degree that will lead 
to the formation of a disk. If this occurs at early stages of the post-AGB phase, it can lead 
to the formation of accretion disk that launches two jets, and by that shapes the PNs. We 
note that the low mass companion, M 2 ~ 0.1M*, will not accrete much mass by itself. 

To summarize this entire section, we did not demonstrate that a backflow forms in the 
PN phase; this will require intensive 3D numerical simulations, because the dynamical range 
from perturbations to a backflow to the center is very large. We only showed here that it is 
quite plausible that dense blobs will flow from the nebula back to the central star, and will 
cause the star outer layers to rotate fast; it might even form an accretion disk. These effects 
will result in an axisymmetrical mass ejection if a VLTP take place. 



3. THE INFLUENCE OF ACCRETION ON A VERY LATE THERMAL 

PULSE 

Various objects have their appearance (in particular, chemical abundances) explained 
as due to a final thermal pulse occuring during the transition from the AGB to the PN 
phase. To name a few examples (see e.g. recent reviews by Werner & Herwing 2006 and 
Kimeswenger et al. 2008): FG Sagittae (the prototypical born-again AGB object), Abell 30, 
Abell 78, CK Vul, PG 1159 (and the class of pulsating white dwarfs named after it), and 
WR-like central stars of PNe. If the above description of PN fall-back is correct, it is worth 
checking if a connection can be drawn between the fall-back and VLTP. 

In general, three types of "late" thermal pulse scenarios can be distinguished (e.g., 
Herwig et al. 1999, Blocker 2003). (1) AGB Final Thermal Pulse (AFTP) takes place at 
the very end of the AGB phase, right before the star leaves this evolutionary stage. (2) Late 
Thermal Pulse (LTP) occurs when a star evolves off the AGB at phase > 0.85 of the thermal 
pulse cycle and it is caught by a TP on the horizontal part of its post-AGB evolution in the 
H-R diagram. (3) Very Late Thermal Pulse occurs still later, when the hydrogen burning 
shell is extinguished and a star is descending along the white dwarf cooling track towards 
lower luminosities. In the VLTP case, the H-rich envelope is penetrated by the pulse-driven 
convective zone which leads to an additional short-lived H-flash. The abundance patterns 
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produced in these three cases are somewhat different, but they all produce H-deficient ejecta 
and H-deficient post-AGB objects. 

Can the accretion of back-falling PN matter, as described in the previous section, affect 
when and how a "late" TP occurs? By definition, it cannot affect an AFTP (no PN yet at 
this stage), and interfering with an LTP seems precluded by the requirement that the fast 
wind should practically cease for fall-back to become possible (Sect. 12.31) . That does not 
happen before moving below the high temperature knee on the HR diagram (e.g., Perinotto 
et al. 2004). Unlike for gravity driven backflows considered by Soker (2001, 2008), this late 
PN backflow cannot prolong the post-AGB lifetime at the horizontal part of the track and 
increase the chances for an LTP (as proposed in Zijlstra et al. 2001 for CK Vul). Note, 
however, that many main sequence late O dwarfs seem to have winds about two orders of 
magnitude weaker than expected from radiative winds theory (see Marcolino et al. 2009 
for a recent update on the topic). If a similar phenomenon would exist among post-AGB 
stars, it would make fall back to the central star more likely, be it before or after the hight 
temperature knee in HR diagram. We will now focus on the VLTP. 

According to sections 12.21 and | 2.3[ fall-back would add 10~ 5 — 10~ 3 M Q of H-rich matter 
to the envelope of a PN central star before a VLTP occurs. Burning of this material can 
be either stable on unstable, depending on the exact conditions (accretion rate, metallicity, 
central star mass and intrinsic luminosity, see Shen & Bildsten 2007, esp. their Fig. 3). The 
relatively high intrinsic luminosity of the young pre- WD (a few hundred L & ) has a stabilising 
effect on H-shell burning and should be taken into account in the stability considerations, 
unlike in classical novae. In the O.6M core mass model used by Perinotto et al. (2004), at 
the time when fall-back is expected to occur, i.e, at the age of 8000 — 10 OOOyrs, the intrinsic 
luminosity is ~ 10 36 erg s -1 . According to the results of Shen & Bildsten (2007), this leads 
to unstable H-shell burning only in a narrow range of accretion rates, 0.3 — 6 x 1O _8 M yr _1 . 

In the case of a stable H-shell burning, the effect of the accreted matter would mainly 
be to prolong the residual H-shell burning and keep the star longer in the evolutionary stage 
not far below the high-temperature knee on the H-R diagram. Accreting 10~ 5 - 10" 3 M Q of 
H-rich matter on a O.6M pre- WD at a stable-burning rate of ~ 10~ 9 M Q (Shen & Bildsten 
2007) would take 10 4 — 10 6 yrs. This is comparable to the duration of a TP cycle (~ 10 5 yrs 
for a 0.6M Q core mass) and processes a similar amount of matter (an AGB star processes 
10 -2 — 1O _3 M between two consecutive TPs, with the higher values characterizing lower 
mass cores, < 0.6M Q ). Effective fall-back would therefore increase the number of post-AGB 
objects experiencing a VLTP above the standard estimate of ~ 10% (Iben 1984). The 
increase would be more significant for the more massive pre- WD objects, as they require less 
mass processed and accumulated in the He-shell for a TP to occur. Notice also that after 
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accretion the star will rotate due to the angular momentum transferred, which may affect 
the flash intrinsic geometry (in addition to constrainig the ejecta, as described in Sect. 12.41) 
and possibly other characteristics. This is different from "normal" late TPs, in which the 
core is rather not expected to rotate rapidly. 

On the other hand, if the fall-back accretion process is slow enough (but not as slow 
as to be stabilized by the residual core luminosity), H-shell burning becomes unstable. In 
this case, fall-back would also switch on a thermonuclear instability on the pre- WD, but in 
the form of H-shell flashes. They would be similar to weak nova eruptions, with the H-shell 
igniting rather than the He-shell (as is in a VLTP). However, from the external appearance 
- new ejecta within an old PN - they could easily be taken for a VLTP. In the past, Iben 
& MacDonald (1986) introduced a term self-induced nova for a final flash of H-burning they 
found to happen late in the evolution of a cooling WD with a very thin He layer buffering 
between the CO core and the H-rich envelope. In their case the ignition was due to diffusive 
mixing of C and H through this thin buffer and occured very late in the WD cooling stage 
(2 x 10 7 yrs after the PN stage, when the luminosity of their model dropped to 0.3L©). The 
fall-back induced flash described here would also be a self-induced nova, but of a distinct kind, 
less violent and much earlier in the evolution, when the pre- WD is still hot and luminous 
and its ejected PN is still around. 



4. SUMMARY 

We have presented a mechanism for a backflow of matter to occur relatively late during 
PN evolution. The proposed fall-back consists of dense blobs that are formed from insta- 
bilities at the discontinuity between the slow and the fast winds of the central star. They 
detach from the dense shell of the slow AGB wind at the beginning of the PN ionization 
stage and fall back into the hot bubble with a velocity of several km s _1 . They can reach 
the PN central star only when the fast wind mass loss rate drops to < 1O _11 M yr _1 . 

The accreted mass is likely to posses enough angular momentum to spin-up the central 
star, and even to form an accretion disk. The estimated accretion rates of M acc ~ 10~ 9 — 
10~ 6 M Q yr _1 , lasting for hundreds to thousands of years, may lead to the formation of an 
accretion disk and to launching of jets, even in the case of single central stars of PNe. The 
accreted mass, ~ 10~ 5 — 1O~ 3 M , would be enough to increase the chance of a VLTP or 
even induce an H-shell flash, leading to a mild self-induced nova within a PN. The VLTP (or 
H-flash) ejecta would exhibit bipolar morphology, offering a binary-independent explanation 
of point-symmetric inner nebular structures prevalent among purported VLTP objects. 
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The estimates provided in this paper do not constitute a proof of the existence of such 
PN fall-back, but they show that it is a reasonable possibility. Numerical simulations of 3D 
radiative hydrodynamics, detailed enough to resolve the blobs, will be required for a definite 
statement. In any case, this scenario, together with other types of post-AGB backflow (as 
described by Soker 2001 and 2008), serves to extend the theoretical argument for the reality 
and importance of backflow in the PN phase. 

This research was supported by the Asher Space Research Institute in the Technion and 
by the Israel Science foundation. 
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